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Abstract 

We discribe a simple way to derive spin correlation functions in 2D Ising model at 
critical temperature but with nonzero magnetic field at the boundary. Local magnetization 



> 



(i.e. one-point function) is computed explicitly for half-plane and disk geometries. 



m 

On 

^ | As is known, the continuous limit of 2D Ising model at critical temperature T = T c 

is described by c = 1/2 conformal field theory which admits a Lagrangian formulation as 
a theory of free massless Majorana Fermi field (ip, ip) with the action 

A=±- f d 2 z[i(>d^+$d z $}. (1) 
X ■ 2tv J 

Here (z, z) = (x + iy, x — iy) are complex coordinates of euclidean plane and d z = 
dxdy. The factor — 1/2tt in this action corresponds to standard normalization of the Fermi 
field, i.e. 

^(z)^(z') = + ... ; = -^—j + - • (2) 

z — z' z — z' 

The order parameter field a(w, w) is related to the above Fermi field non-locally (it creates 

square-root branch cut for ip, ifj) so that the the components i/j(z),i(j(z) satisfy operator 

product expansions 



ijj(z)a(w, w) = (z — w) 2 *y~)(z — w) n a- n cr( 



n=0 
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if)(z)a(w,w) = (z — w) 2 5^(2 — w) n a- n a(w, w), (3) 



n=0 



where a_ n and a_ n are mode operators associated with ip and The field a^a is known 
as "disorder parameter" field \i\ in fact we have 



ao<j(w, w) = —=n(w, w); ao<j(w, w) = —=n(w, w), (4) 
V 2 v 2 



where 



u = e 4 ; u; = e 4 , (5) 
and we have assumed the standard normalization of the fields a and //, i.e. 

a(w, w)a(0, 0) = (ww)~^ x / + tD) / u(0, 0) = (tytD) - ^ x / + (6) 

where / is the identity operator. In fact, the first two subleading terms in (3) are expressed 
in terms of derivatives of //, namely 

a-Kr(w,w) = -^=(4:d w n(w,w)); a_ 2 (r(w,w) = -^=(-dlij(w,w)). (7) 

a_i<r(u;,u)) = -^=(4<9„/i(w, w)); a^(w,w) = ^=(-dl/i(w,w)). (8) 
These equations are easily obtained with the use of explicit expressions 

T(z) = -i(v>0,VO(*); T(z) = -\{^){z) (9) 

for the components T = T zz and T = T 2S of stress-energy tensor in terms of Fermi fields 
ip, Similar equations hold for the products ip/j,, ifjfx; 



ip(z)fi(w, w) = —j=(z - w) 2 [a(w, w) + 4(z - w)d w a(w, w) + -(z — mj) 2 <9 2 (7(10, w) + ...]; 
\/2 3 

to 1 8 

ip(z)fj,(w,w) = —=(z-w)~z [a(w, w)+4(z-w)du,a(w, w) + -(z-w) 2 d%a{w, w) + ...]. (10) 
V 2 3 

This form of operator product expansions makes it possible to derive linear differential 
equations which determine multipoint correlation functions of the fields a and fx(see e.g. 

[I])- 

Conformal field theory of critical Ising model in the presence of boundary (along with 
more general conformal field theories) is studied in Refs[2-4]. It is shown in [3] that there 
are two essentially different possibilities to choose "critical" boundary conditions, the ones 
which preserve conformal invariance. If the boundary B is given by a parametric curve 



B: z = Z(t); z = Z(t), (11) 
where t is some real parameter along the boundary, these two possibilities are 



and 



where 



[e^ip-e^ip] B = (12) 



[e3^ + e*V] B = 0, (13) 



e(t) = jZ{t\ e(t) = j t Z{t) (14) 

are the components of the vector (e, e) tangent to the boundary; in what follows we assume 
that the parametrization (11) of the curve B is chosen in such a way that this vector has 
unit length, i.e. 

e(t)e(t) = 1; (15) 

we also assume the choice of orientation of this curve in which the unit vector (ze, — ie) 
normal to the boundary points inside the domain T> surrounded by £>; B = dT>. From 
a microscopic point of view, (12) and (13) correspond to "free" and "fixed" boundary 
conditions, respectively. Namely, (12) appears when the microscopic boundary spins are 
not restricted, while (13) corresponds to the situation when the boundary spins are fixed 
to be all in the same position, +1 or —1 (so that (13) represents, in fact, two boundary 
conditions, "fixed to +1' and "fixed to —1"). In both cases, due to the conformal invariance, 
one can reduce the case of a generic (one-component) curve B to the "standard" geometry 
where V is the upper half-plane y > 0. Then (12) and (13) simplify as 

[ip — ifj]z=z = for "free" case, (12a) 



[ip + V'] z=z = for "fixed!' case (13a) 

These equations show that the field ip(z) coincides (up to a sign which depends on whether 
we choose "free" or "fixed" boundary condition) with the analytic continuation of ip(z) to 
the lower half-plane. Therefore it is still possible to use (7,8,10) to derive complete set of 
linear differential equations which determine correlation functions of the fields a and fi in 
the presence of boundary [2] . 

The above "free" and "fixed" boundary conditions are "critical" , i.e. they correspond 
to "fixed points" of the renormalization group flow of "boundary interactions", and as- 
sociated field theories with these boundary conditions are conformally invariant. In this 
paper we consider more general boundary conditions which break conformal symmetry. 
Specifically, we analyze the case when boundary spins interact with constant external 
magnetic field Hb- Clearly, this boundary condition "interpolates" between the "free" 
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case at H B = and "fixed to ±1" cases at H B — > ±00, and so it corresponds to "boundary 
flow" from "free" boundary condition down to "fixed" one [5]. This "flow" can be inter- 
preted as "free" boundary condition "perturbed" by "boundary spin operator" a B (t) ■ The 
latter was identified in [3,4] (see also [6]); it can be expressed in terms of Fermi fields as 

a B (t) = ia(t)[e*ij; + e*$](t), (16) 

where a(t) is auxiliary "boundary" Fermionic degree of freedom with the two-point function 
{a(t)a(t')) f ree = \ sign(t — t') [6]. With this, the full action describing the continuous limit 
of T = T c Ising model with boundary magnetic field takes the form 



Ah = 7r- I d 2 z[ipd- z ip + ipd z ip}+ [ {-^-ipip + ]-aa}dt + ih [ a(t)(e^ + e^)(t)dt, (17) 
^ Jv Jb 4yr 2 Jb 

where a = ^a. The first two terms here give the action for Ising model field theory 
with " free" boundary condition (10) and the last term is the "perturbation" describing 
the interaction with the boundary field; it contains dimensional constant h ~ [length] -1 
(appropriately rescaled external field H B ) and breaks conformal symmetry. Note that the 
full action is quadratic, so it is still a free field theory. The fact that the boundary magnetic 
field does not destroy free-field structure of Ising model is very well known (see e.g. [7]). 
Nevertheless, we consider it interesting to observe that although conformal symmetry is 
broken the correlation functions still satisfy linear differential equations. 

The boundary condition for the fields ip, V> in the presence of the boundary field h is 
derived directly from (17), 

(j t +i\)m = (j t -imt), (is) 

where 

m = e^{t)^{Z{t)); m = e*(t)$(Z(t)), (19) 

and 

X = 4nh 2 . (20) 

We want to show that in the presence of boundary magnetic field h the correlation 
functions of a and \i still satisfy linear differential equations. We will demonstrate this 
explicitely for the one-point function (a(z, z))h ("local magnetization") in the simpest case 
of half-plane geometry. So, we assume that the boundary B is the real axis z = z and T> 
is the upper half-plane y > 0. In this case the boundary condition (19) reduces to 

[(d 2 + i\)1>(z) - {d- z - i\)$(z)] z=z = 0. (21) 
This form of the boundary condition makes it explicit that the fields 

X (z) = (d z + i\mz)- x(z) = (d- z - i\)$(z), (22) 
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enjoy the desired property that x(z) coincides with the analytic continuation of x( z ) to 
the lower half-plane. So, for instance, the correlation function (x{z)fj,(w, w)) is an analytic 
function of z on the full z-plane with two square- root branch points at z = w and z = w. 
Taking into account (10) and the asymptotic behavior 

x(z) ~ z~ x as z — > oo (23) 

one can write 

/ / s / >, , si/ si A(w,w) A{w,w) „. , , 

z — w z — w 

Now, one can use the operator product expansions 

X(z)fj,(w,w) = -^(z-w)~%(~<t(w,w) + (z - w)(2d w + i\)a(w, w) + 

+(z - wf{Ad 2 w + 4i\d w )a(w, w) + ...); (25) 

and 

x(z)fJ,(w, w) = -7= - w)~%(—-a(w,w) + (z — w)(2dw - i\)a(w, w) + 
V 2 2 

+ (z - w) 2 (4d| - AiXd^w, w) + ...), (26) 

which follow from (10), to express the coefficients A, A and B in (24) in terms of the 
function (cr(w,w)) and its derivatives. For instance, expanding (24) in powers of z — w 
with the use of (25) we get 

A(w,w) = ^(w-w)*(--){<r(w,w)) h \ (27a) 

s A(w,w) UJ , 1 1 .. , ,. /^„,s 

S («;,«; + V _ = -j= (w - w)H2d w + i\ -)(a(w,w)) h ; 276 

w — w v 2 4w- w 

• 4(t "'" ) - ' i :(tl ,-»,)i(4aJ + (« + ^-)a„ + 



(lW — Vj) 2 y/2 W — W' 

1 iA 11.,... . . 

+ 7J -+77w ^2 )Hw,iv )) h . 27c 

2 u; — w 16 (w — w) z 

On the other hand, in view of (21), the expansion of this function in powers of z — w is 
controlled by the operator product expansion (26). Therefore we have also 

A(w,w) = ^(w-w)^(--)(a(w,w)) h -, (28a) 

s A(w,w) u> , ,1. . ,11...,. . 

B(w,w)+ 2 ' = —(w- w )*(2d iB -i\--r- ){<r(w,w)) h ; (286) 

w — w 4w — w 



A(w,w) U! . .1 ,.„o , . ., 1 .„ 

- =(«;-«;) ^ (4^ + (-4zA+- )^- 



— to) 2 1 w — W 

1 iX 11.,.., . . 

2w — w 16 (to — u;) 2 

Compatibility of (27a,b) and (28a,b) requires that 

(d w + d iB ){<T(w,w)) h = 0. (29) 

This equation just expresses the translational symmetry of the system: the one-point 
function (a(w,w))h does not depend on x = ^(w + w). With this, (28c) becomes identical 
to (27c) and reduces to the ordinary second-order differential equation 



for the function 



a(2Xy) = (a(iy,-iy)) h . (31) 



This can be reduced to degenerate hypergeometric equation and we find that the only 
solution which does not grow exponentially as Y — > +oo is 

a(Y) = C x Y§#(1/2,1,Y), (32) 

where 

1 r°° 

V(a,c,x) = / e-^-^l + t) - - 1 ^. (33) 

r(«) Jo 

The Y-independant constant C can be fixed by comparing the Y — > asymptotic of the 
solution (32) 

a(Y) = Cx -Lyl[_iog(y) + 0(1)] (34) 

'7T 



with the first nontrivial perturbative contribution 

r+L 

(a(iy, -iy))h = h J (a(iy, -iy)a B (x)) free dx + 0(h 3 ) = 

= 2ih(2y)i(\og(2L) - log(2y)) + 0(h 3 ); (35) 

where we introduced an infrared cut-off L and assumed that the field a is normalized as 
in (6). Finally, 

(a(iy, -iy)) h = A*2* (2y)i*(l/2, 1, 2Xy). (36) 
Note that the Y — > oo asymptotic 
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(<7(iy,-iy)) ft ~2*(2y)-* (37) 

of this function coincides with (a(iy,—iy))fi xe d(see [4]), as it should, since h — > oo limit 
corresponds to the "fixed to +1" boundary condition. 

The above analysis can be repeated in the case when V is a disk of radius R, i.e. the 
boundary B is the circle 

Z{t) = ^exp(^); Z(t) = J Rexp(-^). (38) 
Obviously, in this case the one-point function (cr(z, z))h depends only on the combination 



(a{z,z)) h = s(zz/R 2 ), (39) 
and the function s(X) can be shown to satisfy 

rj d 2 ,4(1 + o) 1 , d ,9 1 l-o , 

where o = Ai?. This equation is solved in terms of hypergeometric functions; the solution 
which is regular everywhere inside the disk V is 

Hz, z)) h = gh *W±d (fl/4)-i (1 - X)lF(l/2, 1/2 + o, 1 + o, X), (41) 

where 

X = ^ 9 = ^h 2 R, (42) 

and F(a, 6, c, X)) is the hypergeometric function. Again, we used the first nontrivial order 
of perturbation theory in h to fix the overall normalization factor in (41). In the limit 
when R goes to infinity while the distance from (z, z) to B is kept finite (i.e. X — > 1 with 
Y = XR(1 - X) fixed), (41) degenerates to (32). 

Similar arguments can be used to derive differential equations for multipoint functions. 
However, these equations are more complicated and we do not describe them here. 
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